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ABSTRACT 

This  note  describes  a  cellular  array  algorithm  for  perform¬ 
ing  a  general  class  of  geometrical  operations  on  a  nxn  digital 
image  in  0(n)  time. 
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1.  Introduction 


Cellular  arrays  are  arrays  of  processing  elements  arranged 
in  a  Cartesian  grid,  with  each  processor  ("cell")  connected  to 
its  neighbors  in  the  grid.  Such  an  array,  with  its  cells  opera¬ 
ting  synchronously,  is  a  very  appropriate  multiprocessor  archi¬ 
tecture  for  parallel  image  processing  [1] .  Several  large  cellu¬ 
lar  arrays  have  actually  been  built;  the  largest  to  date  is  NASA's 
Massively  Parallel  Processor  (MPP),  which  is  of  size  128x128. 

Efficient  cellular  array  algorithms  for  most  of  the  basic 
classes  of  image  processing  operations  are  presented  in  [2,3]. 

These  operations  include  local  operations,  where  the  value  of  an 
output  pixel  depends  only  on  the  values  of  the  corresponding  in¬ 
put  pixel  (in  the  same  position)  and  a  set  of  its  neighbors;  and 
statistics  computations,  where  the  value  of  an  output  pixel  is 
(e.g.)  a  weighted  sum  of  all  the  pixel  values  in  the  input  image. 
(The  values  of  a  transform  of  the  image,  such  as  its  Fourier 
transform,  are  all  weighted  sums  of  this  type. ) 

The  purpose  of  this  note  is  to  outline  a  cellular  array  algo¬ 
rithm  for  performing  a  general  class  of  geometric  operations 
("warpings" )  on  an  image.  A  geometric  operation  is  defined  by  a 
coordinate  transformation  of  the  form  x=f(x',y'),  y=g(x',y'), 
where  x,y  and  x' ,y'  are  coordinates  in  the  input  and  output  images, 
respectively.  The  value  of  the  output  pixel  at  (x',y')  depends 
only  on  the  values  of  the  input  pixels  in  a  neighborhood  of  (x,y), 
where  the  relationship  between  (x,y)  and  (x',y')  is  defined  by 
the  pair  of  functions  f  and  g. 


An  MPP  algorithm  for  rotating  an  image  through  an  arbitrary 
angle  is  described  in  [2] .  This  note  generalizes  that  algori¬ 
thm  and  shows  that  it  applies  to  a  large  class  of  geometric 
operations,  including  all  the  standard  operations  composed  of 
translations,  rotations,  reflections,  magnifications,  perspective 
transformations,  etc.  (For  a  discussion  of  image  resampling 
(and  rescaling)  by  a  cellular  array  see  [4].) 


I 


Cellular  geometric  operations 

In  an  arbitrary  geometric  operation,  the  value  of  each 
output  pixel  depends  only  on  the  values  of  a  bounded  set  of 
input  pixels,  but  these  input  pixels  are  at  an  arbitrary 
location  relative  to  the  output  pixel,  and  the  relative  loca¬ 
tion  may  be  very  different  for  different  pixels.  In  order 
to  compute  the  value  of  an  output  pixel ,  we  must  bring  the 
needed  input  values  into  the  output  position;  but  it  is  not 
obvious  how  to  do  this  simultaneously  for  all  the  output  pixels, 
since  the  needed  information  must  travel  along  a  different  route 
in  each  case,  and  this  may  lead  to  communication  conflicts.  Note 
that  for  simple  geometric  operations  such  as  translation,  it  is 
easy  to  shift  all  the  needed  values  synchronously,  without  con¬ 
flict;  but  it  is  not  clear  how  to  do  this  for  an  arbitrary 
operation.  A  brute-force  solution  would  be  to  shift  the  entire 
input  image  in  such  a  way  that  every  value  passes  through  every 

position,  thus  allowing  each  output  pixel  to  collect  the  values 

2 

that  it  needs;  but  this  shifting  would  require  0 (n  )  time  for 
an  nxn  image.  In  this  section  we  show  how  to  do  the  necessary 
shifting  in  0(n)  time  for  a  large  class  of  geometric  operations. 

To  illustrate  the  approach,  suppose  first  that  the  desired 
geometric  operation  has  the  following  property:  Let  denote 
the  set  of  input  pixels  needed  to  compute  the  value  of  output 
pixel  (i,j);  let  denote  the  i^  row  of  the  image,  and  Cj  the 
j  column,  where  l-ri,j5n.  Let  us  assume  that  for  all  i,j, 


the  set 


(U 

i 


)r\Ri 


is  of  bounded  size,  i.e.,  its  size  does  not  grow  with  n.  Note 

that  U  S. .  is  the  set  of  input  pixels  needed  to  compute  the  out- 
i  ID 

x. 

put  values  for  the  entire  j  column;  this  set  itself  is  not  of 
bounded  size,  but  we  assume  that  its  intersection  with  any  given 
row  is  of  bounded  size.  If  this  assumption  is  true,  we  can  use 
the  following  procedure  (compare  [2])  to  compute  the  values  of 
all  the  output  pixels  in  time  0(n).  The  procedure  begins  with 
the  input  image  stored  in  the  cellular  array,  one  pixel  per  cell. 

a)  Cyclically  shift  the  input  image  horizontally.  During 

this  process,  each  cell  (i,j)  "sees"  every  input  value 

in  R^,  and  can  collect  all  those  values  belonging  to  the 

set  (U  S.  .  )HR.  .  Since  this  set  is  of  bounded  size,  the 
i  ID  1 

cell  has  room  to  store  it. 

b)  Now  cyclically  shift  the  stored  values  vertically.  Dur- 

fcll 

ing  this  process,  each  cell  in  the  (arbitrary)  j  column 

"sees"  all  the  values  in  L)[(U  S..)0R.  ]=U  S..  -  i.e.,  it 

1  i  J-D  1  i  ID 

sees  all  the  input  values  needed  to  compute  all  the  output 
values  in  its  column.  It  can  thus  collect  all  the  data 
it  needs  to  compute  its  own  output  value. 

Thus  at  the  end  of  the  process,  each  cell  can  compute  its  output 
value.  The  cyclic  horizontal  and  vertical  shifting  processes 
each  take  time  0(n),  so  that  the  entire  procedure  takes  time  0(n). 
[Cyclic  shifting  is  not  essential  to  the  algorithm;  the  same 


effect  can  be  accomplished,  in  at  most  twice  the  time,  by 
shifting  a  copy  of  the  input  data  first  in  one  horizontal 
direction,  then  in  the  other,  and  similarly  for  the  vertical.] 


The  procedure  just  described  assume  that  (U  S.  .)flR.  is 

i  ij  i 

of  bounded  size;  but  this  is  not  true  for  all  geometric 

operations.  For  example,  consider  the  operation  that  simply 

transposes  the  image;  then  to  compute  each  column  of  output 

values  we  need  a  row  of  input  values,  so  that U  S. .=R. .  In 

i  ID  J 

fact,  consider  the  operation  of  rotation  by  90°  about  the  center 

(i  ,j~);  then  to  compute  the  values  in  C .  ,  we  need  all  the 
u  u  ^0 
values  in  R.  ,  i.e.,  U  S. .  =R.  . 

10  *  ^0  10 

To  see  how  to  handle  these  cases,  note  that  we  can  trans¬ 
pose  an  image  by  simply  shifting  the  input  data  (cyclically)  in 
the  diagonal  direction.  [The  cellular  array  need  not  have  dia¬ 
gonal  neighbor  connections;  to  achieve  "diagonal"  shifting  we 
simply  alternate  horizontal  and  vertical  shifts.]  This  actually 
brings  each  input  value  past  the  output  cell  that  needs  it, 
which  can  thus  store  it,  so  that  when  the  shifting  is  complete, 
all  the  needed  output  values  are  in  their  proper  places.  (Dia¬ 
gonal  shifting  brings  the  input  value  (i,j)  into  the  successive 
positions,  e.g.,  (i+l,j-l),  (i+2 , j-2 ) , . . . ,  modulo  n;  these  are 

just  the  positions  (h,k)  for  which  h+k=i+ j ,  so  that  they  include 
the  position  (j,i).) 

The  use  of  diagonal  shifting  to  transpose  an  image  suggests 
a  generalization  of  our  geometric  operation  procedure,  based  on 


(cyclic)  shifts  in  two  directions,  not  necessarily  horizontal 
and  vertical.  Let  0  be  any  direction,  e.g.,  horizontal,  ver¬ 
tical,  diagonal,  or  more  generally,  a  direction  defined  by  a 
given  periodic  sequence  of  moves  -  e.g. ,  in  the  first  octant 
(0595j),  sequences  of  k  horizontal  moves  separated  by  single 
diagonal  moves.  Let  Rq  denote  a  "0-row",  i.e.,  a  succession 
of  pixels  lying  along  direction  8.  Our  generalized  procedure 
depends  on  the  following  general  assumption:  There  exists  a 
pair  of  distinct  directions  a, 3  such  that  the  set  (  U  S. . )D 

R_  is  of  bounded  size,  for  all  a-rows  R  and  all  3-rows  R„ . 

P  a  p 

If  this  is  true,  we  proceed  by  cyclically  shifting  the  input 
image  along  the  p -rows,  and  allowing  each  cell  to  store  all  the 
values  that  will  be  needed  to  compute  the  outputs  in  its  a -row. 

We  then  cyclically  shift  the  stored  values  along  the  a-rows,  so 
that  each  cell  can  collect  the  data  needed  to  compute  its  own 
output  value. 

The  generalized  assumption  is  evidently  valid  for  all  of  the 
standard  types  of  geometric  operations.  For  translation  and  re¬ 
scaling  we  can  use  the  horizontal  and  vertical  directions;  for 
reflection  and  rotation,  we  use  a  pair  of  directions  that  do  not 
reflect  or  rotate  into  one  another  (e.g.,  for  transposition  we 
use  a  diagonal  direction;  for  90°  rotation,  we  use  a  pair  of  non¬ 
perpendicular  directions).  Perspective  transformations  are  also 
easily  handled,  since  such  a  transformation  can  be  thought  of  as 
a  scale  change  (by  a  different  amount)  along  each  line  parallel 
to  a  given  direction. 


3.  Concluding  remarks 


It  is  known  [1,3]  that  a  cellular  array  can  perform  local 
operations  on  an  n*n  image  in  0 (const)  time,  and  can  perform 
statistics  computations,  and  compute  discrete  transforms,  in 
0(n)  time.  In  this  note,  generalizing  a  result  in  [2] ,  we  have 
shown  that  a  cellular  array  can  also  perform  a  wide  class  of  geo¬ 
metric  operations  on  an  image,  including  all  the  standard  types 
of  operations,  in  0(n)  time. 
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